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ABSTRACT 


This  study  investigates  Magnus  effects  on  a  non-finned  missile  model  which  had  three 
axially  rotating  sections  that  spun  at  different  rates.  Five  spin  cases  are  examined;  the  mid 
section  spinning  only,  die  front  and  rear  sections  spinning  only;  all  sections  spinning  at  die  same 
speed;  all  sections  spinning  in  the  same  direction  with  the  mid  section  spinning  at  an  additional 
10  rev/sec;  and  all  sections  spinning  with  die  mid  section  spinning  in  the  opposite  direction  at  an 
additional  10  rev/sec.  These  five  cases  are  tested  at  three  different  spin  rates  for  two  different 
wind  tunnel  velocities.  Since  no  research  on  differential  spinning  was  found,  potential  flow 
theory  of  a  single  spinning  cylinder  is  expanded  to  the  differential  spinning  cases.  Test  results 
show  that  the  measured  Magnus  forces  and  moments  are  much  smaller  than  die  potential  flow 
predictions  though  most  were  in  die  predicted  direction.  Results  also  show  that  normal  and  axial 
forces  are  not  affected  by  spin  which  agrees  with  potential  flow  theory. 


A  WIND  TUNNEL  STUDY  OF  MAGNUS  EFFECTS  ON  A  DIFFERENTIALLY 


ROTATING  MISSILE 

I.  INTRODUCTION 


In  circulation  theory,  a  spinning  cylinder  generates  a  side  force  in  the  direction 
normal  to  the  free  stream  flow  which  is  commonly  called  a  Magnus  force  (Figure  1-1) 
(7:388).  When  studying  missile  motion,  this  side  force  along  with  die  moment  it  produces, 
can  if  not  taken  into  account,  cause  the  missile  to  drift  off  it's  trajectory.  Magnus  effects 
have  been  studied  by  many,  but  only  in  the  realm  of  a  single  spinning  cylinder.  When  a 
missile  with  three  axially  spinning  sections  rotate  at  different  rates,  the  magnitude  and 
direction  of  Magnus  effects  are  unknown. 


Figure  1-1.  Magnus  Force 


This  experiment  was  based  on  die  Folding  Fin  Aircraft  Rocket  (FFAR)  program 
whose  payload  section  spins  axially  at  a  different  rate  than  the  rest  of  the  rotating  missile. 
The  missile  is  77  inches  long  with  a  2.75  inch  diameter  and  has  a  maximum  rotation  rate  of 
26  revolutions  per  second  (rps)  with  the  payload  section  differential  rate  in  die  range  of  +/- 
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1 S  rps.  When  the  FFAR  program  was  estimating  the  Magnus  effects  on  ifs  missile,  the 
data  used  was  for  a  2D  single  spinning  cylinder  (Figure  1*2)  (6:10)  which  was  not  directly 
applicable  to  this  missile.  The  objective  of  this  experiment  was  to  provide  the  FFAR 
program  with  realistic  Magnus  effect  data  based  on  an  approximate  1/2  scale  model  of  the 
FFAR  as  well  as  get  a  basic  understanding  of  differential  Magnus  effects. 


Velocity  ratio.  — 

LL 

Figure  1-2.  Data  on  a  Rotating  Cylinder  (6:26) 
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EL  TEST  MODEL  AND  INSTRUMENTATION 


Test  Model 


The  test  model  shown  in  Figure  2-1,  was  a  0.47  scale  model  of  the  FFAR.  It 
consisted  of  a  solid  aluminum  shaft  with  three  composit  cylindrical  shells  which  rotated 
axially  at  different  rates.  The  three  shells  reflect  die  nose  cone,  the  payload  section  and  the 
motor  case  section  of  the  FFAR.  In  this  study,  they  are  referred  to  as  the  front,  mid,  and 
rear  sections,  respectively. 

Shells.  The  shells  were  tubular  cylinders  that  were  constructed  by  spiral  wrapping 
fiberglass  and  graphite  sheets.  A  fiberglass  sheet  of  0.09  inch  thickness  was  wrapped  at  30 
degrees  off  axis  then  a  graphite  sheet  of  die  same  thickness  was  wrapped  1 80  degrees 
opposite  to  the  fiberglass.  The  thickness  of  each  shell  was  0.18  inches  which  gave  an 
outside  diameter  of  1.305  inches  and  an  inner  diameter  of  1.125  inches.  The  mid  section 
was  a  length  of  12.375  inches  as  scaled  to  die  payload  section  of  die  FFAR.  The  rear 
section  was  a  length  of  20.0  inches  of  which  0. 5  inch  was  a  brass  cylinder.  The  front 
section  was  a  sharp  cone  shape  of  length  2. 1 5  inches  with  an  indued  angle  of  33.6  degrees 
phis  1 . 1  inch  of  cylinder  length  (Figure  2-1).  This  allowed  the  front  section  to  be  driven 
directly  (discussed  under  Motors). 

Aluminum  Shaft.  The  aluminum  shaft  was  made  from  7075-T7  aircraft  grade 
aluminum  with  a  length  of  28.062  inches  (Figure  2-3).  The  front  of  the  shaft  was 
hollowed  out  to  a  depth  of  about  two  inches  for  die  front  motor  to  sit  in.  An  end  cap  of 


2-1 


2-2 


0.0625  inch  thickness  was  placed  over  the  motor  to  keep  it  from  coming  out  of  the  pocket. 
To  drive  the  other  sections,  one  motor  was  used  in  the  mid  section  and  two  in  the  rear 
section.  For  each  motor,  a  circular  cutout  of  0.3543  inch  radius  depth  and  1.625  inches 
long,  as  shown  in  Figure  2-3.  In  the  mid  section  area  between  the  enclosed  front  motor 
and  die  mid  section  motor,  the  shaft  was  turned  down  to  a  0.4375  inch  diameter.  This  was 
done  to  lighten  the  mid  section  and  move  the  center  of  gravity  toward  the  rear.  To  give 
strength  back  to  this  section,  four  seven  inch  long  balsa  wood  struts  were  attached  in  this 
region.  In  die  area  between  the  two  rear  section  motors,  the  shaft  was  cut  into  an  I-beam 
configuration  to  lighten  the  model  weight  also.  The  last  cut  out  in  the  shaft  was  for  the 
balance.  The  rear  of  the  shaft  was  hollowed  out  to  a  depth  of  4.308  inches  with  a  diameter 
of  0.687  inch.  Inside  of  this  area,  a  balance  receiver  was  placed.  It  was  made  of  brass 
with  a  length  of  2.813  inches,  width  of 0.687  inch  and  a  thickness  of  0.0935  inch.  This 
was  done  because  the  balance  gives  better  readings  when  fitted  against  brass. 
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Motors.  The  motors  were  18  mm  in  diameter,  ironless  rotor  DC  motors  (4:77). 
Attached  to  each  motor  was  a  0.687S  inch  diameter  aluminum  disk  with  a  thickness  of 
0. 156  inch.  A  groove  was  made  in  each  disk  to  allow  a  0.5  inch  diameter  plastic  O-ring  to 
sit  securely  on  the  disk.  The  O-ring  nibbed  the  interior  of  the  shell  and  die  shell  turned 
due  to  the  friction.  To  adjust  the  amount  of  O-ring  contact,  a  set  screw  was  drilled  through 
the  shaft  opposite  each  motor.  This  screw  was  then  used  to  adjust  the  motor  up  or  down 
for  the  O-ring  to  have  optimum  contact  with  the  shell  Since  the  front  section  used  direct 
drive,  O-  rings  were  not  used.  Instead,  a  larger  aluminum  disk  of  1. 125  inch  diameter  and 
0.1875  inch  width  was  used.  It  was  attached  to  the  motor  shaft  by  two  set  screws 
tightened  down  against  it  It  was  then  attached  to  the  front  shell  by  two  other  screws. 

Bearings.  For  each  section  to  rotate  independently,  a  total  of  five  bearings  were 
used;  one  for  the  front  section,  and  two  each  for  the  mid  and  rear  sections  (Figure  2-3). 
The  bearings  were  radial  phenolic  retainer  open  bearings  with  a  width  of  0. 1 56  inch. 

Their  outer  diameter  was  1.125  inches  with  an  inner  diameter  of  0.875  inch  (10:24).  The 
three  furthest  rear  bearings  were  placed  on  the  shaft  and  the  shaft  notched  so  that  the 
bearing  would  not  move.  The  front  two  bearings  were  removable  to  allow  for  model 
assembly. 

Sensors.  To  measure  spin  rate,  internal  reflective  optical  sensors  were  used  in  the 
mid  and  rear  sections.  The  sensors  consisted  of  a  transmitter  and  receiver  in  a  single 
housing  (12:18).  The  mid  section  sensor  was  placed  0.3125  inch  behind  the  mid  section 
bearing  while  the  rear  section  sensor  was  placed  0.094  inch  in  front  of  the  rear  section 
bearing  (Figure  2-3).  To  trigger  the  sensor,  a  white  line  was  painted  in  the  interiors  of  the 
shells  so  that  the  sensor  counted  each  revolution.  The  front  section  could  not  house  a 
sensor  but  was  instead  synchronized  with  the  rear  section. 
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Figure  2-3.  Aluminum  Shaft  Schematic 


Instrumentation 

Force  Balance.  The  balance  used  was  a  0.5  inch  Mark  II.  six  component  strain- 
gage  balance  designed  by  Able  Corp.  (Figure  2-4)  (17:3).  Normal  force  and  pitching 
moment  were  measured  with  two  normal  gages,  N1  and  N2,  which  were  spaced  2. 1  inches 
apart  The  maximum  total  normal  force  allowed  was  1 6  Ibf  and  a  16.8  in-lbf  for  pitching 
moment.  Yaw  force  and  yaw  moment  were  measured  with  the  side  gages,  Y1  and  Y2 
which  were  1.70  inches  apart.  The  maximum  total  yaw  force  allowed  was  10  Ibf  with  a 
8.5  in-lbf  yaw  moment  maximum.  For  axial  force,  the  gage,  AX,  was  allowed  to  measure 
a  maximum  of  5  Ibf.  The  last  gage,  roll  moment,  RM,  had  a  maximum  of  2  in-lbf.  The 
accuracy  of  each  gage  was  "+/-  0.25%  of  maximum  load  or  +/-  0.5%  of  applied  load  when 
compared  with  the  best  straight  line  fit"(17:4).  This  meant  normal  force  accuracy  was  +/- 
0.04  Ibf,  yaw  force  accuracy  was  +/_  0.025  Ibf,  and  axial  force  accuracy  was  +/-  0.0125 
Ibf. 


2-5 


Power  Supplies.  Three  power  supplies  were  used  in  the  experiment.  Two  were 
used  to  power  die  four  motors  inside  the  model.  A  Hewlett-Packard  620SB  Dual  DC 
Power  Supply  was  used  to  run  die  nose  and  rear  sections,  while  a  Hewlett-Packard  6236B 
Triple  Output  Power  Supply  was  used  to  run  the  mid  section.  A  Hewlett-Packard  6205C 
Dual  DC  Power  Supply  powered  the  optical  sensors  which  measured  the  rotation  rates  of 
the  mid  and  rear  sections. 

Multimeters.  Three  multimeters  were  used  to  measure/track  the  voltage  needed  to 
maintain  a  certain  spin  rate  for  each  spinning  section.  The  multimeters  for  the  nose  and 
rear  sections  woe  both  Hewlett-Packard  3466A  Digital  multimeters  while  a  Hewlett- 
Packard  3438A  Digital  multimeter  was  used  for  the  mid  section. 

Counters/Strobe.  Counters  were  used  to  set  die  rotation  rates  of  the  mid  and  rear 
sections.  A  Recal-Dana  1992  Nanosecond  Universal  Counter  was  used  for  the  mid 
section,  while  a  Hewlett-Packard  S316A  Universal  Counter  was  used  for  die  rear  section. 
A  strobe  was  used  to  set  die  front  section  to  the  same  spin  rate  as  the  rear  section. 
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Facilities 


AFTT  5  Ft  Wind  Tunnel.  The  AFIT  5  Ft  wind  tunnel  was  built  in  1919  at  the  old 
McCook  Field  in  Dayton,  Ohio.  It  is  a  wooden  circular  tunnel  that  is  an  open  circuit, 
continuous  flow  tunnel.  It  has  a  contraction  ratio  of  3.7  to  1  and  a  5  ft  closed  test  section. 
The  top  speed  in  the  test  section  is  200  mph  and  is  induced  by  two  counter-rotating  12  ft 
diameter  fans. 

HP  38S2B  Acquisition  Unit.  This  digital  data  acquisition  unit  took  all  readings 
from  the  test  section  except  the  spin  rates  of  the  sections  and  the  atmospheric  pressure 
which  were  operator  entered  into  die  computer.  In  this  experiment,  the  following  voltage 
readings  were  taken:  the  six  balance  gages  (i.e.  Nl,  N2,  Yl,  Y2,  AX,  and  RM),  angle  of 
attack  (alpha),  yaw  angle  (beta),  dynamic  pressure  (Q),  free  stream  temperature,  and  the 
base  pressure  of  die  model. 

Computer  Acquisition  Software.  An  integrated  program  runs  all  aspects  of  data 
collection,  reduction  and  presentation.  Table  2-1  shows  the  program  menu  and  a 
description  of  each  choice.  It  contains  a  library  of  subroutines  which  collect  and  reduce 
calibration,  tare,  and  run  data.  Collection  routines  can  be  modified  to  allow  for  used 
defined  inputs  or  programs  such  as  collecting  spin  rates  in  this  experiment  Reduced  data 
can  be  presented  in  many  forms,  i.e.,  balance  forces,  body  forces,  body  axis,  stability  axis, 
and  wind  axis.  However,  the  reduction  program  was  designed  for  standard  aircraft 
coefficient  reduction,  which  was  a  limitation  for  this  experiment. 
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Table  2-1. 


Computer  Program  Overview 


Setup 

System  Configuration 
Model 

. Writes  the  system  configuration  file 

...Writes  model  parameters  to  a  file 

Tunnel 

...Writes  tunnel  parameters  to  a  file 

Corrections 

...Writes  correction  selections  to  a  file 

Boundary 

..Writes  variable  limits  to  a  file 

System  Calibration 

Gage 

...Reduces  rawdatato  the  calibration  matrix 

Pressure 

...Writes  voltages  from  the  pressure  transducers 

Angle  to  Volts 

Bend 

...Writes  voltages  from  the  angle  transducers 
...Writes  bend  of  sting  vs.  force 

System  Tares 

Average 

Reduce  to  Slopes 
Define  Acquisition 

Run  Acquisition 

..Averages  tare  rawdata  files 

...Reduces  rawdatato  tare  slopes 

...Writes  run.def  files  which  controls  acquisition 

...Selects  and  runs  acquisition  programs 

Run  Reduction 

...Selects  and  runs  reduction  programs 

Pickout 

...Selects  and  presents  reduced  data 

Graph  run 

...Selects  reduced  data  for  graphing 
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m.  THEORY 


When  a  spinning  missile  is  at  some  angle  of  attack  relative  to  the  free  stream  flow,  a 
side  force,  known  as  a  Magnus  force,  is  generated .  This  force  is  named  after  the  physicist, 
G.  Magnus,  who  in  1852,  first  studied  and  measured  this  effect  (8:84).  Most  research  on 
Magnus  effects  assume  that  when  a  circular  cylinder  with  a  spin  rate,  p  (rev/sec),  is  placed 
perpendicular  in  a  uniform  flow  stream,  V,  a  lift  force  develops  (5:2,  13:345)  (Figure  3-1): 

L  =  pVF  (per  unit  length)  (3-1) 

where: 

F  =  circulation 
p  =  density  of  air 

Put  in  kinematic  terms,  the  lift  force  is  the  cross  product  of  the  circulation  and  the  velocity 
(13:346). 
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Eq  (3-1)  is  the  Kutta-Joukowski  Theorem  which  states 


"if  there  is  a  circulation  T  around  an  infinite  cylinder  of  arbitrary  cross  section  placed 
in  a  uniform  stream  V,  the  cylinder  will  experience  a  lift  force  L  =  p  V  F  per  unit 
length  (7:358)." 


This  circulation,  F,  is  defined  as  "the  line  integral  of  the  velocity  around  any  closed  curve" 
symbolically  shown  as: 


V  «dr 


(3-2) 


where 


V  -  velocity  vector 

r  =  vector  radius  from  any  fixed  origin 
c  =  any  closed  curve 

When  circulation  is  applied  to  a  cylinder  such  as  in  Figure  3-2,  it  becomes  (16:267-271): 


r  =J>r  d<9=  2wpr2Jd0 

0  0 

(3-3) 

r  =  4/r2pr2 

(3-4) 

where:  p  =  spin  rate  (rev/sec) 
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Using  Eq  (3-4)  in  Eq  (3-1)  yields  a  lift  force  of: 

L  =  p  V  (4  /r  2pr 2 )  (per  unit  length)  (3-5) 

However,  for  a  spinning  missile,  Eq  (3-5)  is  not  appropriate.  A  spinning  missile's  Magnus 
force  should  be  modeled  by  a  rotating  circular  cylinder  facing  into  the  free  stream  velocity 
not  perpendicular  to  if  s  line  of  symmetry. 

Magnus  Effects  on  a  Single  Cylinder 

Magnus  Force.  A  cylindrical  missile  has  a  coordinate  system  defined  as  in  Figure 
3-3.  The  X-axis  is  positive  out  the  nose,  the  Y-axis  is  positive  out  the  right  side  looking 
from  the  rear,  and  the  Z-axis  is  positive  toward  the  ground.  The  free  stream  velocity  is 
defined  in  the  x-axis  of  the  inertial  system.  The  spin  rate,  p,  of  the  cylindrical  missile  is 
defined  as  rotating  clockwise  as  seen  from  the  re  v  *  me  cylinder.  From  the  definition  of 
circulation  in  Figure  3-2,  the  circulation  of  the  cylindrical  missile  would  be  defined  as  -T 
with  the  vector  along  the  positive  X-axis.  In  this  configuration,  any  force  generated  due 
to  spinning  is  defined  as  a  Magnus  force. 
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Generalizing  to  only  a  cylinder  as  in  Figure  3-4,  at  zero  angle  of  attack,  the  free 
stream  velocity  is  along  the  negative  X-axis  of  the  cylinders  axis  frame.  In  this  situation, 
no  force  can  be  generated  since  the  velocity  and  circulation  vectors  are  both  in  the  same 
axis  and  their  cross  product  is  zero. 


Therefore,  the  only  time  a  Magnus  force  is  generated  is  when  the  cylinder  is  at  some  angle 
of  attack,  a.  When  the  cylinder  is  at  an  angle  of  attack,  the  velocity  vector  has  two 
components  in  the  cylinder's  body  frame,  i.e.  V  sin  a  and  V  cos  a.  It  is  the  normal 
velocity  component,  V  sin  a,  that  is  perpendicular  to  die  circulation  vector  and  can 
generate  a  Magnus  force,  Yp  (Figure  3-5).  The  Magnus  force  generated  is  along  die 
negative  Y-axis. 
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To  define  the  Magnus  force,  Eq  (3-5)  is  translated  to  the  cylinder  system  by  using 
the  normal  velocity,  V  sin  a,  instead  of  the  free  stream  velocity  V,  and  the  circulation  is 
defined  as  -T  instead  of  T.  This  leads  to  the  definition  of  the  Magnus  force,  Yp,  for  a 
single  spinning  cylinder  as: 

Yp--pVsinar  (3-6) 

then  substituting  Eq  (3-4)  into  Eq  (3-6)  leads  to: 

Yp  =  p  V  sin  a  (-4  n  2pr2  )  (per  unit  length)  (3-7) 

For  a  cylinder  of  length,  /: 

Yp  =  />Vsina(-4ar^pr^)  /  (3-8) 

If  dynamic  pressure  is  defined  as  q  =  \p  V2,  the  longitudinal  cross  sectional  area  of  the 

cylinder  is  defined  as  5/  =  2r /,  velocity  ratio  is  defined  as  VR  =  (~— ),  and  for  small 

angles  of  attack,  sin  a  =  cl,  then  Eq  (3-8)  can  be  rewritten  as: 

Yp  =  -2nqSt  VR??  (3-9) 

To  put  it  in  aeronautical  coefficient  form, 


Cv  = 


=  -2n 


(3-10) 


q  St  VR  a 

Magnus  Moment.  The  Magnus  moment  is  the  Magnus  force  multiplied  by  the 
distance  from  the  eg  to  the  Magnus  force  center  of  pressure  which  is  assumed  to  be  in 
front  of  the  eg.  (13:346).  This  distance  is  denoted  by  x.  Taking  the  cross  product  of  the 
distance,  x,  which  is  along  the  positive  X-axis,  and  the  Magnus  force,  -Yp,  the  Magnus 
moment  is  along  the  negative  Z-axis  or: 

np=Y,x  (3-11) 


The  coefficient  form  of  the  Magnus  moment  is  the  same  as  in  Eq  (3-10),  but  with  the 
addition  of  the  cylinder  length. 


3-5 


=  -2  nx 


(3-12) 


C.  = - * - 

"»  qSt  a  l  VR 

Magnus  Effects  with  Partially  Spinning  Sections 

Magnus  Forces.  For  a  cylinder  which  has  two  sections  that  rotate  at  the  same 
speed  as  the  single  cylinder,  each  section's  Magnus  force  would  be  the  same  as  in  Eq  (3-8) 
except  that  the  length,  /,  would  be  the  length  of  the  spinning  section  only  (Figure  3-6) 
denoted  by  ( )f  or  (  ^  for  the  front  or  rear  section,  respectively.  The  Magnus  force  for 
each  section  would  be  a  percentage  of  the  total  length: 

Ypf  =/>Vsina(-4*V  /f)  =  (y)Y„  (3-13) 

Ypr  =/jVsina(-4>r 2pr2  /,)  =  (^)  Yp  (3-14) 


Figure  3-6.  Rear  Section  Spinning  Only 
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However,  their  coefficients  would  be  based  on  the  total  length  or: 


CY  - 


Y. 


Pf 


ri  q St  a  VR 


(3-15) 


Cv  = 


(3-16) 


pr  qS,  a  VR 

Each  coefficient  for  a  partially  spinning  cylinder,  would  be  a  percentage  of  a  single  cylinder 
based  on  the  same  total  length,  l,  or: 


=  (— 

Tpf  y  l '  qSt  a  VR 


Cvf  =Ci)  —  = -2«(*) 


rL 

l 


/  / 

Cy  =(-x) - 1 - =  -2rc(—) 

"  VqS,  crVR  / 


(3-17) 


(3-18) 


Magnus  Moments.  For  the  Magnus  moments,  each  would  be  based  on  half  their 
respective  lengths,  assuming  that  the  eg.  of  the  total  cylinder  is  at  the  point  where  the  two 
sections  join  and  that  the  center  of  pressure  for  a  section  is  at  the  mid  point  of  each  section 
or. 


nP,= 


YPf'f  V. 


n 


2.1 


Y  /  Y  /r2 

Pr  r  _  t  * 


(3-19) 


Pr 


2  21 

The  front  section  moment  would  be  in  the  negative  Z  direction  while  die  rear  section 
moment  would  be  in  the  positive  Z  direction .  The  Magnus  moment  coefficients  are 
determined  the  same  as  in  Eq  (3-12)  or: 


(3-20) 


*"  qSt  a  l  VR  V 


(3-21) 
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(3-22) 


n  I2 

- - - =  T[{~) 

qS{  a  I VR  72 


Magnus  Effects  with  Differentially  Spinning  Sections 

Magnus  Force.  For  a  cylinder  that  has  two  sections  rotating  at  different  rates,  the 
forces  of  each  can  be  calculated  and  then  added  together  to  obtain  the  total  Magnus  force 
(Figure  3-7),  if  it  is  assumed  that  there  is  no  influence  of  the  front  spinning  section  on  the 
rear  spinning  section.  The  Magnus  force  of  the  rear  section  with  a  spin  rate,  p,  would  be 
as  shown  in  Eq(3-14).  For  a  cylinder  with  the  front  section  spinning  in  the  same 
direction  as  the  rear  but  at  an  increased  spin  rate  of  (p+Ap),  the  equation  would  be  similar 
to  Eq  (3-13)  or: 

Ypf  =  p  Vsin  a  [-4;r  2(p  +  Ap)r2/f  ]  (3-23) 
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When  (he  two  forces  are  added  together,  the  resulting  Magnus  force,  is 

equal  to  the  force  for  a  single  spinning  cylinder  of  length  /  plus  an  additional  force 
proportional  to  the  percentage  of  the  front  section  to  the  total  length  for  the  delta  spin  rate: 

=  pV sin  a[-4n  2pr2(/f  +  /,)]  +  pVsin  a(-4n  2  Apr2/,)  (3-24) 

=  Yp  +  AYpf=Yp(l  +  M)  (3.25) 

If  the  front  section  were  spun  in  the  opposite  direction  at  an  increased  spin  rate  of 
-(p  +  Ap),  the  Magnus  force,  YKp+/5p)  would  be: 

Y<,*ap)  =  pVsin  a  [-4*  2pr2(/r  -/, )]  + pVsin  a[-4;r  2(-Ap)r2/f]  (3-26) 

=  Y  -AY  =  Y  (l'-~ I(  - )  (3-27) 

P(Jr'If)  Pf  /  p l} 

In  Figure  3-8,  three  cases  are  shown:  A)  All  sections  spinning  at  the  same  speed;  B)  tire 
mid  section  spinning  faster  and  in  the  same  direction;  C)  the  mid  section  spinning  faster  but 
in  the  opposite  direction. 


A.  B.  C. 


front  =  +p  front  =  (p+A  p)  front  =-(p+Ap) 

rear  =  4p  rear  =  +p  rear  =  +p 

Figure  3-8.  Differentially  Rotating  Sections 
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Note  that  for  case  C,  the  total  Magnus  force  can  be  either  positive  or  negative  depending 
on  the  spin  rate,  p,  and  Ap  chosen.  Each  of  the  forces  would  nondimensionalize  the  same 
as  in  Eq  (3-10)  giving: 


Cy  =-2/r(l+^) 
Y(p+*p)  p  / 

cv  „ 

Ytp^)  /  p/ 


(3-28) 


(3-29) 


Magnus  Moments.  To  calculate  the  Magnus  moments  when  the  front  section  is 
spinning  at  an  increased  spin  rate  of  (p+Ap),  it  is  assumed  that  the  eg.  of  the  cylinder  is  at 
the  meeting  point  of  the  two  sections  and  that  the  center  of  pressure  of  each  section  is  at 
the  mid  point  of  that  section.  The  moment  arm  is  in  the  positive  X  direction  for  the  front 

section  while  for  the  rear  section  it  is  in  the  negative  X  direction.  When  crossed  with  the 
Magnus  force  vector,  the  total  moment  is  in  the  negative  Z  direction.  The  moment,  n^+^j , 

is  equal  to: 


V*)  =  [Ypf  (|  l< )  +  AYpf(K )]  -  Ypr  (j  /, ) 

-lYP(/2  /2iAp/f2) 

~  1  /  ^  f  P  } 


(3-30) 


A  moment  in  the  negative  Z  direction  is  defined  as  a  positive  Magnus  moment.  For  the 
cases  where  the  front  section  is  smaller  than  the  rear  section,  the  moment  generated  would 
become  progressively  less  positive  with  the  possibility  of  becoming  negative. 

When  the  front  section  is  spinning  in  the  opposite  direction  at  -(p+Ap),  the 
moment,  n^p+ap),  would  equal: 


3-10 


(3-31) 


n^,=(-Ypf-AY„)(i/f)-(YPr)(i/r) 
.  Y_ 


f(/r2  +  42 

l  P 


For  cases  where  the  front  section  is  smaller  than  the  rear  section,  the  moment  would  be 
progressively  less  negative  but  not  less  than  die  amount  for  a  single  spinning  cylinder.  To 
obtain  coefficients,  the  Eqs  (3-30  and  3-31)  and  Eq  (3-12)  are  used: 


C 

C 


Model  Predictions. 

Since  this  missile  has  a  pointed  nose  cone  in  addition  to  a  cylinder  body,  the  total 
length  of  the  missile  does  not  take  this  nose  cone  shape  into  account  Some  suggest  that  an 
effective  length  be  used.  Effective  length,  /c gf,  can  be  defined  as  half  the  nose  length, 
plus  die  body  length,  or  (5:4): 

/efif  =  l/2/n  +  *b  (3-  34) 

For  the  model  in  this  experiment,  the  effective  length,  lc fj,  equaled  34.55  inches.  The 
length  of  the  mid  section  was  35.82%  of  the  total  effective  length  and  the  front/rear 
combination  was  64.18%  of  the  total  effective  length. 

Magnus  Forces.  By  using  the  equations  determined  above,  predicted  Magnus 
effects  can  be  determined  for  the  model.  For  the  mid  section  spinning  only,  die  Cypf 
would  equal  (0.3582  Cyp)  or  (-0.7164n)  from  Eq  (3-17).  For  the  front/rear  sections 
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spinning  only,  the  Cypr  would  equal  (0.6418  Cyp)  or  (-1.283671)  from  Eq  (3-18).  The 

for  the  mid  section  spinning  at  a  +Ap,  would  equal:  -2>t(1  +  0.3582  ^- )  from 

Eq(3-28).  When  the  mid  section  is  spinning  in  the  opposite  direction  at  -(p+Ap),  the 

Cy^p+ap)  would  equal:  -2^(0.2836-0.3482— )  from  Eq  (3-29).  Figure  3-9  shows  the 

Magnus  force  versus  angle  of  attack,  a,  graphically,  for  all  five  cases  examined:  mid 
section  spinning  only,  front/rear  sections  spinning  only,  all  sections  spinning  at  the  same 
speed,  the  mid  section  spinning  faster  and  in  the  same  direction,  and  die  mid  section 
spinning  faster  but  in  the  opposite  direction;  based  on  a  spin  rate,  p  =  21  rps,  a  Ap  =  10 
tps.  For  this  spin  rate,  all  Magnus  forces  should  have  a  negative  magnitude. 

Magnus  Moments.  The  Magnus  moment  coefficients  can  also  be  predicted. 

When  the  moment  is  taken  about  the  mid  point  of  the  two  sections  and  the  center  of 
pressures  are  at  the  mid  point  of  each  section,  the  moment  for  both  sections  spinning  in  the 
same  direction  and  at  the  same  spin  rate,  p,  is  zero.  For  the  mid  section  spinning  only, 
C^pf  equals  (-0.128371)  from  Eq  (3-21).  For  the  front/rear  sections  spinning  only, 
would  equals  (0.41 19n)  from  Eq  (3-22).  For  the  case  of  the  mid  section  spinning  a  +10 

rps  faster  in  the  same  direction,  Cn(p+ap)  would  equal  [-  tt(-0.2836  +  0.1283— )]  from 

Eq  (3-32).  When  the  mid  section  is  spinning  at  +10  rps  faster  but  in  the  opposite  direction, 

would  equal  [  0. 5402  +  0. 1283 — )] .  Using  the  same  spin  rates  as  used  in 

Figure  3-9,  Figure  3-10  shows  the  predicted  Magnus  moments. 


3-12 


Figure  3-10.  Predictions  for  Magnus  Moment  Coefficients 
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Single  Cylinder  Research.  Only  limited  research  on  Magnus  effects  has  been 
conducted.  Most  of  the  research  has  been  for  the  infinite  spinning  cylinder  case  as  shown 
in  Figure  3-1  (5:9),  and  no  research  was  found  in  the  area  of  differentially  rotating 
cylinders.  Eq  (3-10),  the  basis  for  the  predictions  in  Figures  3-9  and  3-10,  was  developed 
based  on  2-D  flow  theory,  not  3-D  flow  as  in  the  case  of  a  spinning  missile  at  an  angle  of 
attack.  In  reality,  Magnus  effects  are  caused  by  changes  in  the  boundary  layer  and  vortex 
patterns  and  their  magnitude  and  direction  are  dependent  on;  spin  rate,  angle  of  attack, 
velocity  ratio  and  Reynolds  number.  At  low  angle  of  attack  and  laminar  flow,  the  vortices 
stay  within  the  boundary  layer  which  causes  a  thickening  of  the  boundary  layer  in  the 
direction  of  the  spin  (Figure  3-11).  This  effectively  creates  an  asymmetric  body  shape  that 
generates  a  side  force  in  the  negative  Y  direction.  At  higher  angles  of  attack,  the  laminar 
flow  has  separated  with  two  symmetric  vortices  (Figure  3- 12a).  With  spinning,  the 
separated  flow  shifts  in  the  direction  of  spin  which  also  generates  a  force  in  the  negative  Y 
direction  (Figure  3-12b)  (15:198).  The  magnitude  and  direction  of  the  Magnus  force  is 

also  dependent  on  the  velocity  ratio,  VR  which  is  defined  as  <^~) .  At  low  angles  of 

attack,  the  larger  the  velocity  ratio,  the  larger  the  force  becomes  and  the  closer  the 
transition  region  moves  toward  the  nose  (Figure  3-13)  (5:  20).  At  high  angles  of  attack,  as 
the  velocity  ratio  grows  larger,  the  magnitude  becomes  smaller  and  smaller  until  it  switches 
direction  and  begins  growing  larger  again  in  the  new  direction.  For  this  experiment,  the 
velocity  ratio  was  under  0. 1 5  for  all  test  cases.  And,  though  not  explicitly  stated  in 
Jacobsen,  a  velocity  ratio  of  this  size  and  smaller  is  considered  small.  This  means  that  for 
any  angle  of  attack,  the  Magnus  force  should  be  small  and  in  the  negative  Y  direction. 
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—  Boundary  layer  is  shifted  to 
asymmetric  shape  by  spin. 


- — Cy  (-) 


Figure  3-1 1  Low  Angle  of  Attack  Disturbances  (15:198) 


r~  Separated  region  is  shifted 
„  /  to  asymmetric  shape  by  spin. 

Zt'/ry 

(b)  (ry 

•*—  Cy  (-) 

Figure  3-13.  High  Angle  of  Attack  Disturbances  (15:198) 
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IV.  TEST  SET-UP  AND  PROCEDURES 


In  any  wind  tunnel  experiment,  there  are  preliminaiy  tests  and  checkouts  that  must 
be  done  before  starting  the  test  matrix  in  order  to  insure  that  the  data  obtained  is  accurate 
and  reliable.  In  genera],  the  balance  is  calibrated,  a  tare  run  made  and  checked  out.  Once 
the  preliminaiy  tests  are  completed,  the  test  runs  begin. 

Calibration 


The  force  balance  was  calibrated  to  ensure  that  all  gages  performed  accurately  and 
correctly.  Each  gage  was  calibrated  in  the  positive  and  negative  direction.  Both  are 
necessary  in  order  to  create  the  calibration  matrix  which  was  used  to  eliminate  interaction 
between  the  gages  when  reducing  the  data.  To  calibrate,  the  calibration  body  was  placed 
on  the  balance  with  a  rod  hanging  from  it.  A  weight  was  placed  on  the  rod  and  the 
deviation  from  zero  degrees  was  measured  with  an  inclinometer  in  minutes.  Then  the 
sting/balance  was  raised  back  to  zero  and  the  data  point  taken.  This  procedure  was 
repeated  throughout  each  gage's  range.  The  entire  range  of  each  gage  was  calibrated, 
though  more  data  points  were  taken  at  the  low  end  to  ensure  greater  accuracy  in  that 
region.  From  the  data  points,  the  weight  vs.  voltage  curve  (Figure  4-1),  the  sting  bend  vs. 
angle  of  attack  curve,  and  each  gage's  linear  correlation  were  all  calculated.  The  linear 
correlation  for  each  gage  was  determined  from  die  curves  shown  in  Figure  4-1.  All  gages 
had  excellent  correlations. 
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Figure  4-1.  Gage  Calibration  Curves 


Table  4-1 


Gage  Correlations 


N1 

N2 

Y1 

Y2 

RM 

AX+/- 

Positive 

.99999 

.99997 

.99994 

.99996 

1.0 

.99967 

.99996 

.99998 

.99998 

.99993 

1.0 

Tare  Run 


With  the  model  on  the  balance,  an  angle  of  attack  sweep  from  0-16  degrees  by  2 
degree  increments  was  done  with  no  air  on.  With  this  data,  the  reduction  program  found 
the  weight  of  the  model  and  the  distance  from  the  electrical  center  of  the  balance  to  the  eg. 
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for  each  axis.  To  make  sure  that  this  information  was  correct,  a  test  run  with  no  air  on 
was  done.  This  was  an  a  sweep  from  0-16  degrees  by  2  degree  increments  up  and  4 
degree  increments  down.  The  run  was  reduced  using  the  tare  information.  The  forces  and 
moments  were  evaluated  to  see  how  accurately  the  tare  was.  If  the  tare  was  accurate  and 
the  data  extremely  repeatable,  all  the  forces  would  be  zero  at  every  angle  of  attack.  This, 
however,  was  not  the  case  (and  never  is).  In  order  to  have  the  cleanest  data  possible,  the 
tare  numbers  were  adjusted  until  the  best  solution  was  reached  so  that  the  yaw  force 
would  be  most  accurate.  The  adjusted  weight  was  determined  to  be  1.90  lbf  and  the 
distances  to  the  eg.  were  0.5875  ft  in  the  X  direction,  -0.003843  ft  in  the  Y  direction,  and 
0.035077  ft  in  the  Z  direction  from  the  electrical  center.  Table  4-2  shows  the  error  of  the 
tare.  For  all  angles,  less  than  0.014  lbf  remained  in  the  yaw  force,  less  than  -0.018  lbf  in 
the  normal  force,  and  less  than  -0.056  lbf  in  the  axial  force.  For  normal  and  yaw  forces, 
these  numbers  were  within  theerror  of  the  balance,  i.e.,  +/-0.04  lbf  for  the  normal  force 
and  +/- 0.025  lbf  for  the  yaw  force.  Axial  force  error,  however,  was  more  than  four  times 
the  error  of  the  balance,  i.e.,  +/-0.0125  lbf,  but  this  was  deemed  acceptable  for  this 
experiment. 


Table  4-2. 
Tare  Accuracy 


Alpha 

Beta 

NormFor 

YawFor 

AxialFor  : 

^itchMom 

YawMom 

RollMom 

0.1 

1.07 

-0.0153 

0.0117 

0.001! 

0.0014 

0.0003 

0.0003 

2 

1.1 

-0.0014 

0.0012 

r  -0.01 34| 

0.0017 

0.0001 

0.0001 

3.99 

1.1 

-0.0039 

0.0038 

-0.0146! 

0.0005 

0.0003 

0.0001 

6.1 

1.07 

-0.0177 

0.014 

-0.02751 

0.0022 

0.0011 

0.0001 

8.02 

1.14 

-0.0064 

0.0019 

-0.0331 

0.00281 

00016 

-0.0001 

10.08 

1.14 

0.0096 

0.0026 

-0.042! 

0.0032 

0.0021 

-0.0002 

12.09 

1.1 

0.0071 

-0.00021 

-0.0456] 

0.0028 

0.0017 

0.0002 

14.06 

1.07 

0.0093 

0.0068 

-O.O439! 

0.0014 

0.0024 

0.0003 

16.05 

1.07 

0.0055 

0.0036 

-0.056- 

0.0021 

6.0027 

0.0002 
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Test  Runs 


Preliminary  testing  was  done  to  determine  the  wind  tunnel  speeds  that  caused  the 
least  model  vibration.  Tunnel  speeds  of  114  ft/s  and  162  ft/s  were  decided  on.  By  taking 
the  lowest  spin  rate  which  did  not  vibrate  the  model,  and  the  lowest  tunnel  speed,  a  velocity 
ratio  of  0.042  was  found.  The  rest  of  the  test  matrix  was  determined  by  picking  velocity 
ratios  of  0.07  and  0. 1  and  determining  the  spin  rates  needed  to  obtain  them  (Table  4-3). 

Table  4-3. 

Test  Matrix 


VR  =  0.04 

VR  =  0.07 

VR  =  0.10 

V  =  114  ft/s 

mm 

p  =  33.5  rps 

V=  162ft/s 

mm 

Within  each  of  die  six  above  test  configurations,  five  different  spin  cases  were  done. 
These  were:  a)  mid  section  spinning  only;  b)  front/rear  sections  spinning  only,  c)  all 
sections  spinning  at  the  same  speed;  d)  all  sections  spinning  with  the  mid  section  at  10  rps 
faster  than  the  front/rear  in  the  same  direction;  and  e)  all  sections  spinning  with  the  mid 
section  at  10  rps  faster  than  the  front/rear  in  the  opposite  direction.  From  these  cases,  it 
was  hoped  that  the  contributions  of  each  section,  and  the  influence  of  the  mid  section  on 
the  front/rear  section  could  be  determined. 


4-4 


Test  Set-Up 


The  following  procedure  was  used  when  a  test  run  was  accomplished.  The  flow  of 
data  acquisition  and  the  test  set-up  are  shown  in  Figure  4-2.  1 )  A  barometer  reading  was 
taken  and  manually  input  into  die  Z-248  computer;  2)  A  reading  was  taken  at  zero  angle 
of  attack  and  the  air  off;  3)  The  tunnel  air  was  set  to  the  desired  velocity  and  the  spin  rate 
of  the  scction(s)  set  by  counter  and  strobe  readings;  4)  Once  the  tunnel  reached  the  correct 
velocity,  an  alpha  sweep  from  0-16  degrees  of  attack  at  2  degree  increments  up  and  4 
degree  increments  down  was  begun;  S)  a  data  point  was  taken  at  each  increment  by  the 
HP  3852B  Acquisition  unit  and  the  operator  read  the  spin  rate(s)  off  die  counters  and  input 
them  manually  into  the  Z-248  computer;  6)  after  the  full  alpha  sweep,  the  tunnel  air  was 
turned  off. 
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Once  the  test  matrix  was  complete  and  the  results  reviewed,  eight  test  runs  were 
repeated  in  order  to  test  repeatability  of  the  data.  Repeatability  was  defined  as  good  when 
each  of  the  side  force  gages,  Y1  and  Y2,  were  within  0.3  lbf  for  a  point  within  the  alpha 
sweep.  Four  repeats  were  test  runs  that  did  not  have  good  repeatability  within  the  alpha 
sweep  and  four  were  test  runs  which  did  have  good  repeatability.  Seven  of  die  eight 
repeats  provided  as  good  or  better  data.  Examples  of  good  and  poor  repeatability  for  the 
side  force  gages  are  shown  in  Figures  4-3  and  4-4,  while  Figures  4-5  and  4-6  show  the 
repeatability  of  the  normal  force  gages .  In  Figure  4-3,  each  side  force  gage  had  good 
repeatability,  but  the  total  side  force,  Y1+Y2,  had  variation.  In  Figure  4-6,  the  normal 
gage  repeatability  was  good  even  though  the  side  force  gage  repeatability  for  the  same  test 
run  was  poor  as  seen  in  Figure  4-4. 


0.00  4.00  6.00  12.00  16.00  20.00 

Angle  of  Attack  —  or 


Figure  4-3.  Good  Gage  Repeatability  -  Yaw  Force 
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Mid  Section  Spinning  ot  -4-10  r  ps  Same  0»rection 
p  rr  47.0  rps.  M  — 0.15 
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Figure  4-6.  Poor  Gage  Repeatability  -  Normal  Force 
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V.  DATA  REDUCTION 


Wind  Tunnel  Corrections 


When  doing  experiments  in  a  wind  tunnel,  there  are  a  number  of  corrections  that 
need  to  be  made  to  the  data  in  order  to  take  out  die  effects  that  the  wind  tunnel  itself  may 
impart  to  die  data.  In  this  experiment,  three  were  applied.  They  were:  a  dynamic  pressure 
skew  factor,  a  solid  blockage  correction  to  the  dynamic  pressure,  and  a  buoyancy 
correction  to  the  Axial  coefficient.  A  wake  blockage  correction  to  dynamic  pressure  was 
not  included  nor  was  a  base  pressure  correction  to  the  axial  coefficient.  Not  enough 
information  was  found  to  accurately  calculate  these  factors. 

Skew  factor.  This  factor  accounts  for  the  difference  between  where  the  dynamic 
pressure  is  measured  in  the  tunnel  and  the  dynamic  pressure  in  the  test  section  (11:142). 
The  factor  for  this  tunnel  is  1.019  or: 

Q-^QnK^l-019  (5-1) 

Solid  Blockage.  By  placing  a  model  in  the  tunnel,  the  tunnel  area  where  air  can 
flow  is  reduced.  From  Bernoulli's  equation  and  the  continuity  equation,  die  velocity 
flowing  over  the  model  increases.  This  increase  is  considered  constant  over  the  model  and 
is  called  "solid  blockage".  It  is  calculated  based  on  the  model  volume  and  tunnel  area 
(11:353,364-6)  or: 
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where: 


K( model  volume) 

£*  ~  (5“2> 
=  0.00032 


K  =  the  body  shape  factor  (0.96  for  a  body  of  revolution) 

C  -  the  tunnel  test  section  area  (19.635  sq.  ft  for  this  tunnel) 
This  leads  to  a  new  dynamic  pressure  equal  to: 


q  =  q«**  (1+  **)*  (5-3) 

Buoyancy.  This  factor  accounts  for  the  variation  in  static  pressure  due  to  boundary 
layer  increase  as  the  flow  moves  downstream.  The  boundary  layer  increase  decreases  die 
pressure  along  the  test  section  creating  a  negative  pressure  gradient.  This  pressure  change 
is  in  effect,  making  the  tunnel  section  smaller  downstream  which  causes  the  streamlines  to 
squeeze  together  and  increases  the  drag  on  the  model.  This  drag  increment  is  defined  as: 


AC 


AJ8 


dp  body  volume 


qsr 


) 


(08629) 


(5-4) 


where: 


( — )  -  slope  of  longitudinal  static  pressure  curve 
dl 


This  coefficient  increment  is  subtracted  from  the  axial  force  coefficient  (11:350-3,  362-4). 
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Magnus  Effects 


Magnus  Force.  Though  the  wind  tunnel  acquisition  program  reduced  the  data  to 
force  and  moment  coefficients,  it  used  die  standard  aircraft  nondimensionalization,  q  Sf, 
where  q  is  the  corrected  dynamic  pressure  and  Sr  is  the  cross  sectional  area.  This  was  not 
an  appropriate  way  to  nondimensionalize  the  Magnus  force  coefficient,  since  it  was  also  a 
function  of  the  spin  rate,  p  and  normal  velocity,  V  sin  a.  Instead,  the 
nondimensionalization  of  Eq  (3-10)  was  used: 


)a 


By  using  these  parameters,  the  velocity  ratio,  the  model  length,  and  the  normal  component 
of  the  free  stream  velocity  are  taken  into  account  However,  Eq  (3-10)  is  for  the  case  of  a 
cylinder,  not  a  missile.  To  apply  the  equation  to  this  experiment,  it  must  be  modified  to: 


where: 


Yr« 


)a 


(5-5) 


Ptfr  =  0.6418(^-^)  +  0.3582(pm) 

Y„  =  Yp-N, 

This  equation  now  takes  into  account  the  nose  cone  of  die  missile  and  the  variations  in  spin 
rates  due  to  the  differential  as  well  as  die  differences  in  spin  rates  between  the  front  and 
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rear  sections  at  angles  of  attack.  However,  by  using  the  angle  of  attack,  a,  in  the 
denominator,  the  coefficient  tends  to  fluctuate  greatly  at  small  angles  (a  <  4  degrees),  but 
at  higher  angles,  the  data  tended  to  collapse.  Therefore,  this  variation  at  small  angles  was 
deemed  acceptable. 

The  measured  yaw  force  which  should  have  only  contained  the  Magnus  force  in 
this  experiment,  did,  in  reality,  contain  an  additional  yaw  force  due  to  a  constant  yaw  angle, 
p,  of  1.10  +/-0.05  degrees.  This  yaw  force,  Ng,  was  also  in  the  negative  Y  direction  and  in 
effect  added  to  the  Magnus  force.  If  not  subtracted  out,  the  Magnus  force  appeared  larger 
than  it  actually  was.  From  the  geometry  in  Figure  5-1,  this  force, was  equal  to  (9:966): 

N,  -  N0  sinf  (5-6) 

where: 

Nq  =  the  normal  force  without  spinning 


Figure  5-1.  Geometry  of  Induced  Beta  Yaw  Force  (9:966) 
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This  force  was  calculated  for  each  angle  of  attack  and  subtracted  out  of  the  measured  yaw 
force  to  give  a  corrected  yaw  force,  Yj^.  Table  5-1  shows  this  beta  induced  yaw  force  for 
the  two  wind  tunnel  test  velocities. 


Table  5-1. 

Beta  Induced  Yaw  Forces  (lb#) 


alpha 

Ns-V=120 

CVJ 

* 

z 

2 

■003228 

0.010554 

4 

■0.00463 

-0.00898 

6 

-001697' 

-0.01349 

T1 

■0.01325 

-0.02807 

10 

•0.01237 

l~-0.02665 

12 

-0.01517 

-0.03172 

14 

-0.01958 

-0.03567 

16 

-0.02617 

-0.04095 

From  potential  flow  theory,  it  is  known  that  at  zero  angle  of  attack,  no  yaw  force 
should  be  measured.  However,  due  to  small  model  vibrations,  a  yaw  force  was  measured 
at  zero  angle  of  attack.  This  tended  to  spread  the  data  results.  In  order  to  better  analyze 
the  data,  this  measured  yaw  force  at  zero  angle  of  attack,  was  considered  a  constant  value 
and  subtracted  out  of  each  angle  of  attack  for  each  test  run.  In  essence,  driving  each  test 
run  through  zero.  In  some  cases,  this  force  at  zero  angle  of  attack  was  exceptionally  large 
and  if  subtracted  out,  unduly  shifted  the  entire  test  run.  In  these  cases,  the  test  run  was  not 
zeroed.  These  cases  are  noted  with  an  *  in  die  figures  containing  them.  In  die  Appendix, 
the  measured  forces  as  well  as  the  coefficients  are  provided. 
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Magnus  Moment.  In  order  to  subtract  out  the  beta  induced  yaw  force,  the 


following  equation  was  used: 

%  =  nine_a(1  y)  (5-7) 

where: 

x  -  the  distance  from  the  balance  center  to  the  eg.  of  the  model  in  the  X  direction 
=  0.5875  ft 

y  -  the  distance  from  the  balance  center  to  the  eg.  of  the  model  in  the  Y  direction 
=  0.003843  ft 

A  =  measured  axial  force  (lbf) 

To  obtain  the  Magnus  moment  coefficients,  Eq  (3-12)  was  modified  to: 


n 


p« 


cir 


(5-8) 


In  keeping  with  the  procedure  used  to  calculate  Magnus  forces,  any  measured  moment  at 
zero  angle  of  attack,  was  subtracted  out  of  the  test  run  expec  t  in  the  cases  where  it  was 
exceptionally  large  as  described  for  the  Magnus  forces  earlier. 


Normal  and  Axial  Coefficients 

Normal  and  Axial  force  coefficients  should  not  change  due  to  the  spinning  and 
therefore,  the  standard  nondimensionalization  of  q  ST,  was  appropriate  and  used. 
However,  since  die  model  was  at  a  constant  yaw  angle  when  the  model  was  spinning,  a 
small  Magnus  force  was  generated  in  the  normal  direction.  Since  the  yaw  angle  was  very 
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small  (~  1  degree),  this  additional  force  would  be  very  small  when  compared  to  the  true 
normal  force.  For  this  reason,  it  was  not  subtracted  out  of  the  measured  normal  force. 
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VL  RESULTS  AND  DISCUSSION 


Magnus  Force  Effects 

Mid  Only.  Front/Rear  Only,  and  AB  Sections  Spinning.  In  Figures  6-1,  6-2,  and 
6-3,  there  appears  to  have  been  some  interaction  between  the  mid  and  rear  sections.  In 
Figure  6-3  (all  sections  spinning),  the  magnitude  of  the  coefficient  did  not  appear  to  be  the 
two  individual  cases  (Figures  6-1  and  6-2)  added  together  which  was  what  the  potential 
theory  predicted.  Instead,  all  three  cases  had  about  equal  magnitudes.  This  shows  that 
there  was  interaction  between  the  mid  and  front/rear  sections.  Interaction  means  that  the 
vortices  coming  off  a  section  interfere  with  die  flow  over  the  next  section  which,  in  effect, 
decreases  the  magnitude  of  that  rear  section.  Since  the  predicted  values  were  calculated 
based  on  no  interaction,  this  would  account  for  the  predicted  values  being  larger  than  the 
measured  values. 

AH  Mid  Same  +10  ros.  Mid  Opp.  +10  rps  Sections  Spinning.  When  Figures  6-3, 
6-4,  and  6-5,  were  analyzed,  there  were  definite  differences  in  the  Magnus  force 
coefficients.  When  the  mid  section  was  spinning  at  +10  rps  in  the  same  direction  (Figure 
6-4),  the  Magnus  force  appeared  to  be  the  same  magnitude  as  in  the  all  sections  spinning  at 
the  same  speed  case  (Figure  6-3).  This  could  mean  that  the  interaction  between  the  mid 
and  rear  sections  effectively  canceled  out  the  increased  magnitude  of  the  mid  section  due  to 
the  increased  spin  rate.  Or,  the  increased  Magnus  force  due  to  the  differential  was  too 
small  to  be  measured  accurately.  However,  in  Figure  6-6,  when  the  mid  section  was 
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spinning  in  the  opposite  direction  at  a  +10  rps,  there  was  a  noticeable  change  in  the 
Magnus  force  coefficient  In  theory,  the  magnitude  generated  should  be  small  but  not 
positive,  however,  in  this  case,  the  magnitude  of  the  force  was  positive.  When  the  mid 
section  was  spinning  opposite  to  the  front/rear  sections,  the  mid  section  generated  enough 
force  to  totally  cancel  out  the  front/rear  generated  force  and  still  have  some  magnitude  left. 
What  can't  be  determined  here  was  the  effect  of  the  differential.  To  determine  this,  testing 
of  the  mid  section  spinning  in  the  opposite  direction  at  the  same  spin  rate  as  the  front/rear 
sections  needs  to  be  done.  An  important  result  of  this  case  was  that  this  type  of  spinning 
could  effectively  cancel  Magnus  force  generation  which  would  eliminate  the  problem  of 
trajectory  drift  that  a  Magnus  force  introduces. 
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Figure  6-1.  Magnus  Force  Coefficient  -  Mid  Section  Only 
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Figure  6-2.  Magnus  Force  Coefficient  -  Front/Rear  Sections  Only 
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Figure  6-3.  Magnus  Force  Coefficient  -  All  Sections  Same  Speed 


Figure  6-4.  Magnus  Force  Coefficient  -  Mid  Section  +10  rps  Same  Direction 
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Figure  6-5.  Magnus  Force  Coefficient  -  Mid  Section  + 10  rps  Opposite  Direction 


Magnus  Moment 


Mid  Only  and  Front/Rear  Only  Spinning.  In  Chapter  n,  to  predict  the  Magnus 
moment  it  was  assumed  dial  the  center  of  pressure  for  each  section  was  at  ifs  midpoint.  In 
Figures  6-6  and  6-7,  this  assumption  is  shown  to  be  wrong.  Though  the  moments  were 
small  in  Figures  6-6,  and  6-7,  they  had  about  the  same  positive  magnitude.  In  order  for  the 
moments  to  be  positive  when  their  respective  Magnus  forces  were  negative,  the  model's 
center  of  pressure  must  be  slightly  behind  the  eg  which  was  located  in  the  rear  section  but 
not  at  it's  midpoint.  Also,  die  moment  was  essentially  constant  which  means  that  the  center 
of  pressure  did  not  move  with  increased  angle  of  attack  since  the  Magnus  forces  were  also 
essentially  constant. 

AH  Mid  Same  +10  rps.  Mid  Opp  +10  rps.  From  Chapter  n,  for  all  sections 
spinning  at  the  same  speed,  die  moment  was  taken  about  the  eg  and  was  zero.  In  Figure 
6-8,  for  this  case,  the  trend  for  the  moment  was  a  slightly  positive  magnitude.  This  means 
that  the  center  of  pressure  was  behind  the  eg  as  was  the  case  in  Figures  6-6  and  6-7. 
However,  in  Figure  6-9,  when  the  mid  section  was  spinning  at  +10  rps  faster  in  the  same 
direction,  the  moment  was  approximately  zero.  This  means  that  the  center  of  pressure 
moved  forward  to  the  eg  since  the  forces  for  both  cases  were  essentially  equal  (Figures  6-3 
and  6-4).  Finally,  in  Figure  6-10,  when  the  mid  section  was  spinning  oopositc,  the 
moment  was  slightly  negative.  Since  the  Magnus  force  for  this  case  was  a  positive 
magnitude,  the  center  of  pressure  must  also  be  behind  the  eg.  slightly  as  in  the  previous 
cases.  Overall,  the  moments  were  very  small.  The  differential  spinning  in  either  direction 
did  not  seem  to  greatly  affect  the  moment  generated. 


6-9 


Figure  6-6.  Magnus  Moment  Coefficient  -  Mid  Section  Only 
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Figure  6-7.  Magnus  Moment  Coefficient  -  Front/Rear  Sections  Only 
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Figure  6-8.  Magnus  Moment  Coefficient  -  All  Sections  Same  Speed 
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Figure  6-9.  Magnus  Moment  Coefficient  -  Mid  Section  +10  rps  Same  Direction 
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Figure  6-10.  Magnus  Moment  Coefficient  -  Mid  Section  +10  Opposite  Direction 


Normal  Force 


In  theoiy,  the  normal  force  on  a  spinning  missile  should  not  be  affected  if  it  is  not 
at  a  yaw  angle  because  there  would  be  no  cross  flow  velocity  to  generate  an  additional 
normal  force.  In  this  experiment,  there  was  a  constant  yaw  angle  of  ~1  degree,  but  it  was 
determined  to  be  negligible.  To  determine  if  the  normal  force  had  increased  in  magnitude, 
the  normal  coefficient  of  the  no  spin  case  was  compared  to  coefficients  obtained  from  a 
missile  computer  modeling  program  called  Missile  DATCOM.  This  program  adapts  die 
USAF  Stability  and  Control  DATCOM  Methodology  to  missile  shapes.  It  is  typically  used 
to  predict  aerodynamic  performance  for  preliminary  designs  based  on  data  gathered  from 
previous  missile  programs.  For  this  experiment,  die  physical  characteristics  of  the  model 
were  loaded  into  the  program  and  the  aerodynamic  coefficients  predicted.  It  must  be 
noted,  however,  that  this  model  had  a  length  to  diameter  ratio  of  about  26.S  which  is  out  of 
the  range  of  typical  missiles.  Therefore,  the  DATCOM  predictions  were  only  estimates 
based  on  extrapolations.  In  Figures  6-1 1  through  6-15,  the  DATCOM  predicted  value  was 
consistently  higher  than  the  measured  normal  force.  Though  there  was  some  variation  in 
the  coefficients  which  could  have  been  due  to  die  beta  induced  Magnus  force,  die  curves 
were  all  essentially  the  same.  In  Figure  6-16,  the  normal  coefficients  for  the  no  spin  cases 
were  also  compared  to  the  predictions.  These  verify  that  DATCOM's  prediction  was  too 
large  for  this  model.  However,  die  no  spin  cases  were  the  same  as  the  spinning  cases  and 
could  be  compared  with  die  spinning  cases.  For  example,  when  looking  at  10  degrees 
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angle  of  attack,  the  normal  coefficients  throughout  the  six  cases  (Figures  6-1 1  through  6- 
16)  were  essentially  a  magnitude  of  one.  This  verified  the  fact  that  the  normal  coefficient 
was  not  affected  by  spinning. 


6-15 


6-16 


Figure  6-11.  Normal  Force  Coefficient  -  Mid  Section  Only 


Figure  6-12,  Normal  Force  Coefficient  -  Frort/Rear  Sections  Only 


Figure  6-13.  Normal  Force  Coefficient  -  All  Sections  Same  Speed 


Figure  6-14.  Norma!  Force  Coefficient  -  Mid  Section  +10  rps  Same  Direction 
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Figure  6-15.  Normal  Force  Coefficient  -  Mid  Section  +10  rp°  Opposite  Direction 


Figure  6-16.  No  Spin  Normal  Coefficient 


Axial  Force 


As  in  the  previous  section,  the  Missile  DATCOM  program  was  used  to  generate 
axial  force  coefficients  to  compare  to  the  test  results.  Figure  6-17  shows  the  DATCOM 
results  versus  die  results  for  the  no  spin  cases.  The  measured  axial  force  coefficients  had 
much  greater  variation  than  did  the  normal  force  coefficients.  However,  in  Chapter  IV, 
Tare  discussions,  die  axial  force  error  could  be  as  large  as  0.04  lbf.  As  with  the  normal 
coefficient,  the  predicted  values  were  larger  than  die  actual.  In  Figures  6-18  through  6-22, 
the  axial  force  coefficients  are  shown.  In  all  cases  except  for  the  mid  section  spinning 
opposite  at  +10  rps,  the  trend  was  the  same:  at  zero  angle  of  attack,  the  coefficient  had  a 
magnitude  of  ~0.3  +/-  0. 1  and  then  gradually  decreased  toward  zero  magnitude  at  ~  10 
degrees.  This  is  consistent  with  the  general  trend  of  an  axial  coefficient  for  a  missile  as 
shown  in  Figure  6-23,  taken  from  Chin's  book  on  missile  design  (2:18).  In  Figure  6-22, 
when  the  mid  section  was  spinning  opposite,  die  trend  was  opposite.  Instead  of  decreasing 
with  L  rger  angles  of  attack,  the  magnitude  increased.  This  effect  was  probably  due  to  the 
vortices  of  the  mid  and  rear  sections  coming  off  in  opposite  directions  and  interfering  with 
each  other. 
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Figure  6-17.  Axial  Force  Coefficient  -  No  Spin 
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;e  Coefficient  -  Mid  Section  Only 


00 
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Figure  6-19.  Axial  Force  Coefficient  -  Front/Rear  Section  Only 
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Figure  6-20.  Axial  Force  Coefficient  -  All  Sections  Same  Speed 


Figure  6-21.  Axial  Force  Coefficient  -  Mid  Section  +10  rps  Same  Direction 
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Figure  6-22.  Axial  Force  Coefficient  -  Mid  Section  +10  rps  Opposite  Direction 


o#  9cr 

ANGLE  Of  ATTACK,  a 

Figure  6-23.  Typical  Axial  Coefficients  (2:18) 


vn.  CONCLUSIONS 


Although  definitive  results  were  not  obtained,  general  trends  were  found.  In  most 
cases,  the  potential  theory  and  Missile  DATCOM  predictions  were  too  large.  This  was  due 
to  the  assumptions  made  in  potential  flow  theory  and  extrapolations  made  for  the 
DATCOM  results.  It  was  found  that  there  was  flow  interaction  between  the  mid  and  rear 
sections  which  decreased  the  magnitude  of  the  Magnus  force.  This  also  tended  to  cancel 
out  the  effect  of  the  spin  rate  differential  in  the  same  direction  since  in  Figures  6-3  and  6-4, 
the  magnitudes  were  the  same.  However,  with  a  spin  differential  in  the  opposite  direction, 
the  mid  section  seemed  to  dominate  and  canceled  out  the  Magnus  force  generated  by  the 
rear  section  'Figure  6-5).  Though  the  Magnus  forces  varied,  the  moments  they  generated 
were  quite  small  and  the  center  of  pressures  for  the  five  cases  did  not  vary  much.  The 
center  of  pressure  stayed  slightly  behind  the  eg.  and  did  not  vary  with  angle  of  attack. 
When  the  normal  and  axial  force  coefficients  were  examined,  they  essentially  confirmed 
that  spinning,  whether  differentially  or  not,  did  not  change  their  magnitudes.  Overall,  these 
results  show  that  for  the  Folding  Finned  Aircraft  Rocket  program,  the  Magnus  effects 
should  have  a  very  small,  if  any,  effect  on  the  FFAR's  trajectory.  But  if  possible,  the 
payload  section  should  be  spun  opposite  to  the  front/rear  sections  so  that  no  Magnus 
moment  is  generated. 
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Vm.  RECOMMENDATIONS 


1 .  Though  general  trends  could  be  seen,  the  forces  measured  in  this  experiment  were  not 
very  repeatable.  This  may  have  been  due  to  the  wind  tunnel  velocities  that  were  used.  For 
further  study,  tests  should  be  done  in  a  higher  speed  tunnel.  This  could  eliminate  some  of 
the  noise  in  the  data  as  well  as  better  reflect  the  speeds  the  FFAR  flies  at.  However,  at 
higher  speeds,  model  vibration  might  be  a  greater  problem  so  that  preliminary  testing  at 
various  speed  would  need  to  be  done  first. 

2.  The  front  and  rear  sections  were  supposed  to  spin  together,  but  it  was  not  possible  to 
construct  this  model  to  have  them  internally  connected  to  insure  this.  More  accurate 
results  might  be  measured  if  the  design  was  changed  to  enable  these  two  sections  to  be 
connected. 

3.  To  better  understand  the  effect  of  the  mid  section  spinning  faster  than  the  rest  of  the 
model,  increasing  the  mid  section  spin  by  a  percentage  versus  a  set  rate,  i.e.  +10  rps,  might 
yield  more  collapsible  data. 

4.  The  possibility  of  using  air  to  turn  the  sections  should  continue  to  be  explored.  If 
higher  spin  rates  are  desired,  air  pressure  may  be  the  only  way  to  obtain  them  in  such  a 
small  diameter  model. 
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APPENDIX  :  DATA  TABLES 
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